Exact Groundstates of Rotating Bose Gases close to a Feshbach Resonance 
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We study the groundstates of rotating Bose gases when interactions are affected by a nearby 
Feshbach resonance. We show that exact groundstates at high angular momentum can be found 
analytically for a general and realistic model for the resonant interactions. We identify parameter 
regimes where the exact groundstates are exotic fractional quantum Hall states, the excitations of 
which obey non-abelian exchange statistics. 
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Dilute atomic gases are fascinating systems with which 
to investigate the unusual properties of Bose-condensed 
systems when subjected to rotation. They offer the pos- 
sibility of entering a regime of high rotation rate where 
novel uncondensed groundstates can appear 0. At suf- 
ficiently high angular momentum, the groundstate of a 
system of bosonic atoms interacting by contact repulsion 
is the bosonic Laughlin state[l|,0. In addition, numer- 
ical studies 0j S 13 have indicated that as the angular 
momentum is increased many other strongly-correlated 
groundstates appear before the Laughlin state is reached. 

In order to facilitate the experimental observation of 
the strongly-correlated groundstates, it has been pro- 
posed to use a Feshbach resonance to increase the energy 
scale of two-body interactions^. In addition, the Fesh- 
bach resonance induces an energy dependence to the two- 
body scattering amplitude and leads to a non-zero con- 
centration of bound molecular states of the atoms 8, 9] . It 
is therefore possible that new strongly-correlated ground- 
states may form|^0|- 

In this paper we show that exact groundstates at high 
angular momentum can be found for a general and re- 
alistic model of atoms interacting via a Feshbach res- 
onance. We show that the Laughlin state of atoms 
is the exact groundstate of this model over a range 
of parameter space. Furthermore, we identify parame- 
ter regimes where the exact groundstates are strongly- 
correlated atom/molecule mixtures, which are directly 
related to the unconventional fractional quantum Hall 
states prop osed by Moore and Read|0 and by Read and 
RezayifTij (sometimes referred to as the "Pfafhan" and 
the "parafermion" states). These many-body states have 
very exotic physical properties, including non-abelian ex- 
change statistics of their quasiparticle excitations [lolll]| . 
Exact results are extremely rare in the field of strongly- 
correlated many-body systems, and provide a crucial un- 
derpinning to its study. Our results establish the impor- 
tance of these exotic strongly-correlated states for rotat- 
ing atom gases in a realistic physical situation. 

We consider a system of bosonic atoms and molecules 
confined to a harmonic trap with cylindrical symmetry 
about the z-axis. We denote the natural frequencies of 
the trap by u}\\ and (jJ_l in the axial and transverse di- 
rections. The Hamiltonian we consider is the same as 
that used to explore the effects of a Feshbach resonance 



in other contexts0,0|. We write the Hamiltonian in the 
non-rotating frame in the form 

H ^ Hk + Hp + Hi + {huji_ + fej||/2)iV + huj±L (1) 

defining 
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where dr ^ and rhr^ are the (bosonic) field operators for 
atoms and molecules. 

Hk is the kinetic energy (relative to the zero-point 
energy) in a frame rotating at about the z-axis, with 

K = — \-ihV - Aa{r)f + -MujIz^ 
2M ^ aWJ I 2 II 

where Aa — Mu!±z x r and A„i = 2Mu!±z x r. The 
eigenstates of ha and hm are the products of symmetric- 
gauge Landau level states (in the xj/-coordinates) and 
simple harmonic oscillator states (in z). The single- 
particle states for atoms and molecules differ only in the 
natural quantum lengthscales: in the Landau level states 
the "ma gnetic lengt hs" for the atoms and molecules are 
£a = ^Jh/(2Muj±:) and = £a/V^', in the simple 
harmonic oscillator states, the "oscillator lengths" are 
Wa — y^h/ (Mijj||) and Wm — Wa/V^. (In each case, the 
l/\/2 reduction is due to the molecule/atom mass ratio 
of 2.) 

Hp describes atom-atom scattering via the Feshbach 
resonance|^: g is the strength of hybridisation; the de- 
tuning A is the energy difference between the lowest en- 
ergy state of one molecule and the lowest energy state 
of two atoms in the rotating frame'l2!|, which can be 
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tuned continuously with an applied magnetic field. Hi 
describes elastic atom-molecule and molecule-molecule 
scattering. We shall consider only the case of repulsive 
bare interactions Uaai Uam, Umm ^ and positive detun- 
ing A > 0. However, note that we make no assumptions 
concerning the relative sizes of the interactions and the 
level spacings of the confinement potential. 

The Hamiltonian commutes with both the total num- 
ber of nuclei (in atoms and molecules) 



where 



N : 



\aiar 



2mtmr'\ dPr 



and the total angular momentum (in units of K) 

We shall look for the groundstates of H at fixed N and 
L. The terms in Q proportional to N and L are then 
just constant energy offsets, and we need only consider 
the action of HK^ Hp and Hj. 

The crucial parameter controlling the nature of the 
groundstates of a system of atoms at large angular mo- 
mentum {L '3> N ^ 1) is the filling fraction^. In the 
present situation of an atom/molecule mixture, we define 
the total filling fraction in terms of N and L by 
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Denoting the mean 2D density of nuclei in both atoms 
and molecules by nx = fia + 2n„i, where fia^m are the 
expectation values of the atomic and molecular 2D den- 
sities, then the total filling fraction may be written 

where the individual filling fractions Q for the atoms and 
molecules alone are Pa = "■a(27r£^) and = ^m(27r^m). 
The factor of 4 appearing in this formula is a common 
feature of quantum Hall states containing bound pairs 

la 

We now turn to look for the exact groundstates of the 
Hamiltonian at fixed total N and L. Substantial progress 
can be made by noting that, owing to the contact nature 
of the atom-atom interaction and the atom-molecule hy- 
bridisation, the atomic field operators in Hp can be elim- 
inated in favour of an atom pair operator 

^ a(t)a(t) . 

The resulting quadratic form can be diagonalised to give 
A- 



Hp 



A+ 



Jrhr — aPr 



i (a + ± y(A~CW2)2TV) 



a 



2 

- .g/\/2 

= A_ - A 



Hp is positive semi-dcfinitc provided A_ 
equivalent to the condition 
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< Uaa 



> 0, which is 



(2) 



Under this condition the reduction of energy of a pair of 
atoms arising from hybridisation with the molecular level 
is less than the energy cost from the repulsive contact in- 
teraction to bring the two atoms together: there remains 
a net repulsion between pairs of atoms. 

We shall identify exact groundstates of the Hamilto- 
nian, under the condition Our strategy is to con- 
struct many-particle states (for a given N and L) which 
are simultaneous zero-energy eigenstates of Hk, Hp, and 
Hj. Since each of these terms is positive semi-definite, 
the resulting states are guaranteed to be groundstates of 
the total Hamiltonian for given N and L. 

We start from a purely atomic state |'I')a (for which 
iTir\'^)a = 0), and which is an eigenstate of N and L. 
Consider the (unnormalised) state \'^) formed by 



R 



i?|*>a 

exp J Prihld^r 



(3) 
(4) 



Since R commutes with both N and L, the number of 
nuclei and the total angular momentum of j^") are both 
equal to those of the underlying atomic state \^)a- 

The operator R induces hybridisation of pairs of atoms 
in |^')a into molecules, in such a way that the new state 
© is annihilated by the term containing A-|_ in Hp. This 
follows from noting that [rhrjR] = {a/P)RPr, so 



(3mr - aPr ] i?|^'>a = R (/3to^) |*> 







since rhr\'if)a — 0. We can describe the residual effects 
of Hp on the states (PJ by deriving the action of the 
operator R^HpR within the subspace of purely atomic 
states. Using TOr|5')a = 0, wc find 



RTHpR = A_ 1 



a 

/32 



atdiOo 



(5) 



where O = : exp J P^Pr(Prj : (the colons denote 

normal ordering). We shall refer to the operator jSJ as 
an "effective atomic Hamiltonian". We emphasise that 
it is not a unitary transformation of Hp. However, its 
spectrum does allow one to deduce all eigenstates of Hp 
which are of the form ^ , which is sufficient for our pur- 
poses. The residual interaction is an effective two-body 
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FIG. 1: Processes that generate effective (a) three-body and 
(b) four-body interactions between the atoms. A single (dou- 
ble) line denotes the propagator for an atom (molecule). 

contact repulsion between the underlying atoms (with a 
strength reduced from C/oa/2 by the hybridisation). 
Similarly, the "effective atomic Hamiltonian" from Hi 

is 




+ {fj'^jHtoiarfd^r. (7) 
I 



The atom-molecule and molecule-molecule interactions 
generate effective three- and four-body contact interac- 
tions for the underlying atoms. These arise from the 
processes indicated schematically in Fig. ^ 

We now consider the action of the kinetic energy op- 
erator Hk on the states It is helpful to expand 
the field operators in the eigenstates of ha.m'- we la- 
bel these states by the quantum numbers n, m and 
where n = 0,1,2,... and m = — n, —n -I- 1, . . . are the 
Landau level index and the angular momentum quan- 
tum number (describing motion in the xy-plane), and 
q — 0,1,2,... is the oscillator index of the axial mo- 
tion. We then define ajj^^ = / a\.{r\n,m,q)a, and 
ml^„^g = j <Pr m'l{r\n,m,q)m- The states \n,m,q)m 
and the states \n,m,q)a are each assumed to provide a 
complete basis in which arbitrary wavefunctions can be 
expanded (over all lengthscales for which the contact in- 
teraction Hamiltonian is assumed valid). We may then 
write 

/ PrlTll (fr = ^ Pnmqmi„,^ (8) 

where 



Plmq= i «nimigi42m292 / a("l'"^l'9ll^)'i("2,"l2,92|r')(r|n,TO,g)„ I (9) 



ni,mi,gi n2,ni2,92 



We now insist that all atoms in the state |^')a are in the 
lowest Landau level and lowest axial oscillator state: that 
is, in states with n = q = 0. In this case, the "effective 
atomic Hamiltonian" W HkR vanishes identically. This 
surprising result follows from the fact that the n = q = 
wavefunctions for a molecule, (r|0,m,0)m, which are 
proportional to 

(x + z2/)™e-(^'+'^')/4^™e-^'/2»- 

provide a complete set of states for the product of 
any two atomic wavefunctions with n = q — Q, 
(r|0, mi, 0)a(r|0, m2, 0)a, which is proportional to 

{X + ^y)™i+™2g-2(.^+„^)/<g-2.V2-^ . 

The completeness relation holds since m can span all 
possible values of mi -\- TO2 and since the arguments 
in the exponentials are identical due to tm — ^a/V^ 
and w,n — Wa/V^- Consequently, from (|8I9I) . the ac- 
tion of / PrThl.d'^r on any state in which all atoms have 
n = g = is equivalent to the action of J2m -Pomo'^iomo- 



In the new state ^ molecules are only in the lowest Lan- 
dau level and lowest axial oscillator state {n = q = 0): 
the state (O is therefore annihilated by Hk- We empha- 
sise that choosing to construct |^')a from states in which 
atoms are only in states with n — q ~ does not impose 
further restrictions on the parameters in the model iQ); 
we make this choice because the exact groundstates are 
of this form in the situations we study. 

We are now in a position to deduce exact groundstates 
of the Hamiltonian ^ at particular values of N and L. 
States of the form based on an underlying atomic 
state \'^)a in which all atoms have n ~ q = are zero- 
energy eigenstates of Hk- To identify zero-energy eigen- 
states of Hp and Hj, we further require that the un- 
derlying atomic wavefunction be annihilated by the 
effective two-, three-, and four-body contact interactions 

Zero-energy eigenstates of the two-body contact inter- 
action (and therefore also of the three- and four-body 
interactions), can be constructed in the lowest Landau 
level for any L > N{N — 1) by choosing \^)a to be the 
bosonic Laughlin state (at L = N{N — 1)), or the Laugh- 
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lin state with additional quasiholes (at L > N{N — 1)). 
The LaughUn ket state of N atoms, |^'L)aj may be con- 
structed from its first quantised wavefunction by 

N N 
\^L)a = n Kn*L({n})n«tjO) (10) 

i=l j=l 
*L({n}) = n 

i<j=l 

where rji = {xi + iyi) and |0) is the particle vacuum. (We 
omit normalisation factors.) 

For L = N{N — 1) the Laughlin state is the only state 
in the lowest Landau level that is a zero-energy eigenstate 
of the two-body contact interaction. The state — 
R\'^L)a is therefore the non-degenerate groundstate of 
the model Hamiltonian ^ aX L = N{N— 1), i.e. at total 
filling fraction vt = 1/2. Note, however, that since |^'L)a 
has no amplitude for two atoms to coincide, consequently 
-Pr|*L)a = and = i?|*L)a = |*L)a. The cxact 

groundstate atvT — is the bosonic Laughlin state for 
atoms mUji) . We emphasise that this result is valid for ar- 
bitrarily strong atom-atom, atom- molecule and molecule- 
molecule interactions. The only restrictions we impose 
are that g^/A < Uaa, and A,Uaa,Uam,Ujnm > 0. The 
pairing correlations between atoms suggested by Bhon- 
gale et al.^ do not appear in the groundstate under these 
conditions. 

It is particularly interesting to consider tuning the pa- 
rameters (for example by varying the detuning A) such 
that g^/A = Uaa- Then A_ = and the effective two- 
body repulsion (O vanishes: the effective atomic Hamil- 
tonian consists only of the repulsive three- and four-body 
contact interactions Ht)l7|) . 

Exact zero-energy eigenstates of these interactions can 
be constructed in the lowest Landau level for L > 
N{N/2~l), from the Moore-Read ("Pfafhan") state[U 
The Moore-Read state is the unique zero-energy eigen- 



state of a three-body contact interaction at L = N{N/2 — 

l)i3- Hence, we establish that l^'Mi?) = R\'^MR)a is 
the non-degenerate groundstate of the Hamiltonian 
at L = N{N/2 — 1), which corresponds to z/t = Ip^Slj. 

If Uam/Urnm ^ ff^/(^^w'a^a) ' t^^n the effective three- 
body contact interaction © is negligible compared to 
the effective four-body interaction when vt ~ 0{1). 
It is then of interest to set Uam — and find exact 
groundstates of the four-body contact interaction alone 
ITI). This can be achieved in the lowest Landau level 
for L > N{N/3 - 1) by use of the /c = 3 Read-Rezayi 
( "parafermion" ) wavefunction[Tl] . From the uniqueness 
of this state at L = N{N/3 - 1) [13, we estabhsh that 
|v[/j,^3) = -R|^fc=3)a is the non-degenerate groundstate of 
(HJ at L = iV(iV/3 - 1), corresponding to vt = 3/2. 

Note that the exact groundstates which are formed 
from an underlying Moore-Read state (at i^t — 1) or 
fc = 3 Read-Rezayi state (at vt — 3/2) are strongly- 
correlated atom/molecule mixtures (since Pr\'^MR)a 7^ 
and Pr|^fc=3)a 7^ 0). These atom/molecule mixtures will 
exhibit all the exotic properties that have been predicted 
for the Moore-Read and Read-Rezayi states of atoms 
alone, including non-abelian statistics for their quasipar- 
ticle excitations 0, 0|. Although these states provide 
exact descriptions of the groundstates only at g'^/A = 
Uaa (and Uam = for the Read-Rezayi state), since 
these states have excitation gaps they should continue 
to describe the groundstates at i/y = 1 and i^t = 3/2 
over a range of parameter space. Indeed, these states 
have been found to provide accurate descriptions of the 
groundstates of a purely atomic system (A oo) at 
— l[li^ll| and vt = 3/2 y; it is therefore possible 
that they provide accurate descriptions of the ground- 
states at these filling fractions over the range .g^/A — Uaa 
to g^/A = 0. 
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